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Abstract. We give a new proof of the trace formula for regular graphs. Our ap- 
proach is inspired by path integral approach in quantum mechanics, and calculations 
are mostly combinatorial. 



O ! 1- Introduction 

O 

CN ■ The famous Selberg's trace formula first appeared in [I]. On a compact hyperbolic 

surface it relates the eigenvalue spectrum of Laplace operator to the length spectrum 
of closed geodesies. A version of this formula for finite regular graphs was obtained by 
Ahumada j2] (cf. also Ihara [3]). 

Trace formulae are known to have many implications. For instance, they can be 
considered as nonabelian generalizations of the Poisson summation formula. In case of 
finite graphs, since one can find the eigenvalue spectrum of Laplacian for a given graph 
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{Sj I explicitly, the trace formula lets one find the numbers of closed geodesies of any length 

(see ((HH) ) ■ In physics trace formulae indicate the cases when semi-classical evaluation 
of the path integral for state sum of a quantum free particle in some background is 
exact. Selberg's formula is also known to bear much resemblance to Riemann-Weil 
formula in number theory. 

The original proof of the trace formula for regular graphs (|5U|) was in the framework 
of "discrete harmonic analysis". We propose another way to derive it, inspired by the 
path integral approach in quantum mechanics [Hj. We consider the trace Z&(t) = tr e tA 
as a state sum of the quantum free particle living on the graph. We rewrite it as a 
sum over closed paths, which is a discrete version of the usual path integral over loops 
>• . for quantum mechanical state sum. Then we divide the set of closed paths into classes 

of homotopically equivalent paths. There is a class of contractible paths, and one 
homotopy class for each closed geodesic on the graph. We explicitly calculate the con- 
tribution of each homotopy class to the state sum, thus rewriting it as a contribution of 
contractible paths plus sum over "long" geodesies (of nonzero length) of contributions 
of their individual homotopy classes. This is analogous to the stationary phase calcula- 
tion of the path integral. Geodesies serve as stationary points of the action in the space 
of loops. Homotopy class of a geodesic serves as a neighbourhood of the stationary 
point. Thus we arrive to the known trace formula for regular graph, with a specific, 
physically relevant, choice of test function for eigenvalue spectrum of Laplacian. 
We wish to thank P. Zograf who inspired this work. 
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2. Notations and definitions 

Let T be a finite regular connected non-oriented graph with vertices of valence q+1 > 
2 with no multiple edges and no edges connecting a vertex with itself. Denote by V(r) 
and E(T) the set of vertices of T and the set of edges respectively. Let |T| = #V(r) 
be the number of vertices. Further denote the space of complex-valued functions on 
vertices by Fun(r) = C v( - r \ A basis function (vector) \v > associated with vertex v 

l 



equals 1 on v and on the other vertices. We further adopt the quantum-mechanical 
notations and denote the transposed basis vector by < v\ — \v > T . We call the set of 
vertices connected to v by edges its link and denote it Lk(i>). 

The averaging operator T : Fun(T) — >■ Fun(r) acts as follows: for / G Fun(r) 

(1) (Tf)(v)= f( v ') 

v'eLk(v) 

The Laplace operator A on T is defined by 

(2) (A/)(«)= fW)-Mv)f(v) 

u'eLk(v) 

where val(t> ) is the valence of v. Since we consider a regular graph T, A differs from T 
by a multiple of identity: A = — (q + 1)1 + T where 1 is the identity map 
Fun(T) — > Fun(T). The matrix of the averaging operator is just the adjacency matrix 
of the graph: < v'\T\v >= 1 iff and v' are connected by an edge and otherwise. The 
diagonal elements of T are zero. 

The physically interesting quantity is the trace of the heat kernel (the state sum) 

(3) Z A (t) = tr exp(tA) = e^ 1 )* Z T (t) 
where 

(4) Z T (t) = tr exp(tT) 

in 

It turns out that Zx(t) is more convenient for our calculation than Z\{t). If {A-,}^ 
is the set of eigenvalues of T then 

|r| 

(5) Z T (t) = J2^ t 

j'=i 

If we introduce the spectral density function for T 

|r| 

(6) p(s) = J2^-\ J ) 

j'=i 

we may express Zx(t) as a Laplace transform of p(s): 

/oo 
p(s)e st ds 
-oo 

Eigenvalues Aj are known to satisfy — q — 1 < Xj < q + 1. Moreover, q + 1 is always 
an eigenvalue, while — q — 1 may be an eigenvalue and may be not. If it is, then the 
distribution of eigenvalues is necessarily even p(—s) = p(s). 



3. Sum over paths 

Let us define a closed path of length I as a sequence of vertices (v i, . . . , vi) such that 
for every j = 1, . . . , / the Vj is connected to Vj+i by an edge (we identify vi+i with 
Vi). We will usually omit the word "closed" in the following, since all paths, walks, 
trajectories etc. will be supposed to be closed. Denote the set of paths by P and the 
length of a path p G P by \p\. We also denote the number of closed paths of length / 
by p;. It is convenient to identify paths of length with vertices of Y. 
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Lemma 1. 



This expression may be viewed as a discrete version of path integral over loops for 
the state sum, with j^jj being analogue of the measure e~ s on loops. We give two 
different explanations of (JHJ). The first one is more lengthy, but done in the spirit of 
usual derivation of path integral representation in quantum mechanics. The second is 
absolutely straightforward and evident. 

3.1. First proof of Lemma 1. Let us evaluate Zx{t) in the following manner: 
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(9) Z T (t)= <v\e tT \v>= Jim £ < v \(l + ±T) N 

v£V(r) ^ uev(r) 

= E < ^(l+^HK >< v N \(l+^T)\v N ^ >■■■< V2 \(l+^T)\ Vl > 

vi,...,v N eV(T) 

we are summing here over all sequences of N vertices v\, . . . ,v^. Notice that matrix 
elements < fj+i|(l + jrT)\vi > equal 1 if v i+ i = Vf, 4 if v i+ i and t>j are connected 
by an edge; and otherwise. Let us call a walk of length N a sequence of vertices 
(vi, . . . ,t>7v) such that each pair of successive vertices Vj,Vj + i are either connected by 
an edge or coincide. Let Wn be the set of walks of length N. For a walk w G Wn 
denote the number of values of j for which vj and Vj + \ are connected by an edge by 
\w\. 

The only nonzero terms in the last line of (JHJ) are those with the sequence w = 
(v i, . . . , vn) being a walk. For these terms the summand is (t/N)^. Thus we have 



(10) Z r (f)=Km ^(t/iV) H 

w&Wn 

This is also a sort of discrete path integral representation for Zt- To transform it to 
the form (JHJ), we need a projection ttn '■ Wn — > -P which leaves only those vertices in 
a walk for which Vj ^ Vj+i, and forgets the others. For a walk w G Wjv the result 
of projection n N (w) is a path of length Each path of length I has C l N walks as 
preimages under n N (C l N is a binomial coefficient). So 

(11) Zr(t)= lim J^CW (t/JV) w 

TV— >oo ^— ' 



Using 

( 12 ) lim 77T7 = A 
we obtain (JHJ). □ 

3.2. Second proof of Lemma 1. One can arrive to (JHJ) in a more straightforward 
way: we may just expand the exponent in definition of Zj<{t) in a Taylor series in 
variable t: 

(13) z T (t)=tre tT = J2-trT l 

1=0 



then 

(14) tiT l = <v 1 \V\v l ><v l \V\v l „ 1 >---<v 2 \V\v 1 > 

vi,...,vi<EV(T) 

the terms in this sum with (vi, . . . , Vi) G P equal 1, all the others vanish; hence 

(15) tr T l = Pl 
and we obtain (jHJ). □ 

4. Sum over geodesics 

We use the term "closed trajectory" for equivalence class of closed paths under cyclic 
permutations of vertices along the path. So a trajectory is a path with information on 
the starting point forgotten. An elementary homotopy is a transformation of trajecto- 
ries of the following kind: 

(Vi, . . . , Vj, . . . , V n ) H-> (Vi, . . . , Vj, v', Vj, ...,v n ) 

where v' G Lk(v J ). Two trajectories are called homotopic if they can be connected by a 
chain of elementary homotopies (with arrows either forward or backward). The shortest 
representative in a homotopy class is called a geodesic trajectory (or just geodesic). 
An alternative definition of geodesic trajectory is as a trajectory satisfying Vi ^ Vi +2 
for all i. Denote the set of all geodesics on T by G. Two paths are called homotopic 
if their trajectories are homotopic. If 7 = (v±, . . . , v n ) is a trajectory of length n then 
its r-th power is defined as a trajectory of length rn obtained as 7 walked around r 
times: Y = (vi, . . . , v n , . . . , Vi, . . . , v n ). A trajectory 7 is called primitive if it is not 
a (non-unit) power of any trajectory. A geodesic trajectory with one of its vertices 
chosen as a starting point is a geodesic path. A path homotopic to path of length 
is called contractible. We call geodesics of length short or trivial, and geodesics of 
length > long. 

We proceed now to the calculation of contribution of contractible paths to (JBJ) (one 
may also call it the contribution of short geodesics). 

4.1. Contribution of contractible paths. Let us denote T the covering tree for F 
and call some point C G f the center. The function dist on vertices of the covering 
tree dist : T — > No returns the minimal number of edges one must pass to reach given 
vertex form the center. Any contractible closed path on T can be lifted to a closed path 
on f (and all closed paths there are contractible, since f is a tree) and we adjust the 
lift so that it start and ends in C. Since each edge passed in one direction by a closed 
path on r must by passed in the opposite direction, the length of the path must be 
even. Denote by P^k^) the set of closed paths on T of length 2k starting and ending 
in C; we also need a subset P2k{^) C i"2fc(r) consisting of closed paths not returning 
to C (except the starting point and the end point). 

Recall a concept of Dyck path of length 2k (see e.g. |4J): it is a sequence of integers 
(«i, «2, • • • , «2A;+i) with «i = «2fc+i = 0, ttj > and c^+i = a.; ± 1. We denote the set 
of Dyck paths of length 2k as D 2 k] j^Dik = Cat^ (the k-th Catalan number). There is 
a projection nj? : P2k(X ) ~~ ¥ D 2 k-2- It acts as follows: 

(16) (vi =C,v 2 ,..., v 2 k, v 2k +i = C) i-> (dist(t> 2 ) - 1, ... , dist(w 2 fc) - 1) 

The number of preimages for any Dyck path under irj? equals (q + l)q k ~ 1 since there 
are q + 1 choices to make the step from v± — C to v 2 ; q choices for each step, increasing 
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dist; steps, decreasing dist are done uniquely (since for any vertex v ^ C of T one edge 
from it leads inward, while the q others lead outward). So we have obtained that 

(17) #£»(F) = (g + l^Cat*-! 

The generating function for the numbers of paths #P 2 fc(r) is obtained as a simple 
modification of the usual generating function for Catalan numbers: 

(18) r P {8) = f; #P 2fc (r) ^ = (i + g- 1 ) 1 "^ 1 " 4 ^ 2 

For numbers of paths that may pass through the center we obtain 

(19) 2(l-(g + l)V) 

And hence we obtain the contribution to Z?{t) from contractible paths: 

(20) mou = in • g = in ■ JL rfs e » 1^(1) 

where real part of A is greater than real parts of all singular points of integrand, as 
usual for inverse Laplace transform. The factor of |T| in front is due to the fact that a 
contractible path can start from any vertex of V (we remind that |T| denotes the number 
of vertices in V). Further evaluating the integral we wrap the contour of integration 
around the cut s G [— 2y/q, 2y / g]. Thus we proved 

Lemma 2. The contribution of contractible paths to is 



•s 
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(2D [M«u=iri-ViV se * (o+i) 2 -* 2 

In other words, we obtained a contribution to the spectral density of T on T from 
contractible paths: 



fon\ r ( w in ^ + 1 V 74 ? - s 2 

(22) [p(s)J CO n = |r| 



2tt (g + l) 2 -s 2 
in the interval s E [— 2-Jq, 2-Jq\ and outside it. 

4.2. Contribution of long geodesies. Suppose we have a (long) primitive geodesic 
7 of length I on V. To calculate the contribution of its homotopy class to Z?{t) we 
need to find the number of paths of length k homotopic to 7: p 7i £. Since a path is a 
trajectory with some point on it chosen as a start/end point, p 7i fc = k t 7) fc where t 7j ^ 
is the number of trajectories of length k homotopic to 7. Note that this is not true for 
non-primitive 7. 

We may find the numbers t 7) & using the following combinatorial construction. If 
7 = (vi,...,vi) (with the periodic condition t> = v{), any trajectory that can be 
contracted to 7 can be represented as a closed contractible path from v\ to V\ never 
going along the edge V\ — v ; then a step v 1 — t> 2 ; then a closed contractible path from 
i>2 to i> 2 > never returning along the edge t> 2 — v\, then step ?j 2 — ^3 and so on. All in all 
it is I closed contractible paths with one direction prohibited and I unit steps. Thus 
we find the generating function for t 1 y. 
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(23) ^ traj -\s) = KkS k = s^(F P (s)) hl 

k=h\ 



is 2 



where the superscript (traj.) indicates that we are counting trajectories, jFp{s) is the 
generating function for the numbers of contractible closed paths with one direction 
prohibited: 

(24) M*) = T 1 - 1 ">^ 
and hence 

(25) F\ tra] -' 



^p(s) 2qs 2 



1 - y/1 -Aqs 2 



hi 



2qs 

For the numbers of paths homotopic to 7 we have 

(26) = EW' ^ .|^'W ^ (^Jf^ 

Now we would like to calculate the numbers of paths homotopic to non-primitive 
geodesic 7 = (7') r , that is a primitive geodesic 7' passed r > 2 times. It turns out that 
if we carry out the scheme above in this case, every path becomes calculated r times. 
For any geodesic 7 denote A (7) the length of the primitive geodesic 7 is power of. If 7 
is primitive itself, we set A (7) = |7|. Thus for any geodesic we have 

(07) T"-"-U.) - AiZ) ( L ~V L - A " sl ' 

{ ' 7 2qs 

Now we have all the information to write down the contribution of a long geodesic 
7 to Z T (t): 

A+ioo 



1 pA+ioa 1 1 

(28) [Z r (t)] 7 = — ds e st -F^ aths \- 

2m J A-ioo s s 



1 f A+t0 ° dse st Ml) (s-^^Aq 



2™ J A _ lO0 ^/s 2 - Aq I 2g 



The last integral reduces to the modified Bessel's function of the first kind Ji 7 i. So we 
deduced 

Lemma 3. For every long geodesic 7 G G the contribution of its homotopic class in P 
to d$J) equals 

(29) [Z T (t)] 7 = A(7)g-^ l/2 /| 7 |(2^) 

Collecting together (|2*Tj) and (|2HJ) we obtain the full trace formula: 

Theorem (Trace formula for regular graphs) . Let T be a finite connected regular 
graph of valence q + 1 > 2 with \T\ vertices, without multiple edges and edges connecting 
a vertex to itself; let T be the averaging operator on T and Z T (t) = tr e tT with t a 
complex variable; let G be the set of long closed geodesies on T; for each 7 G G, \^\ is 
the length of '7 and A^) is the length of the underlying primitive geodesic: for 7 = (7')' 
with primitive 7' we set A(7) = |7'|. Then 

(30) Z T (t) = |r| • [ 2V ~ Q ds e* f Aq ~ ** + J2Ml)Q~ hl/2 U^t) 



another useful form of the same result is 



( 31 ) z T{t) = in • »±i * «- ^^p^ + 1 SP , .-npvw 

where gp; is the number of geodesic paths of length I. To pass from (J30|) to (|3*T|) 
one must notice that the number of geodesic paths corresponding to a given geodesic 
trajectory 7 is A (7). 

We may interpret (fHUj) in terms of spectral density p(s) 

Corollary 1. 

(32) p(s) = \F\. q -±l ( _ g2) + ^ A(j) q -^^M= ■ fl(4g - , 2 ) 

2tt (g + l) 2 -s 2 ^ 7Tv/4g-s 2 



where 



(33) 7}(x) = cos(Z arccos(x)) = ~ ((x + Vx 2 - 1)' + (x - Vx 2 - I)') 

is the Chebyshev polynomial of the first kind of degree I. 



The sum on the right of (j32j) is to be understood in the generalized function sense 
(it does not exist in the ordinary sense since numbers gp z grow too fast). Notation 9 
is used for unit step function. A remarkable fact is that although each geodesic gives 
a smooth contribution to p(s) with support on the interval s G [— 2*/g, 2y/q] (with 
singularities on the endpoints), the sum of all contributions is a generalized function 
with support on eigenvalues of T, scattered across a wider interval s G \—q — 1, q+ 1]. 

We may invert (}3*2*j) in a sense to reproduce the numbers of geodesic paths gp ; from 
the spectrum of averaging operator: 

Corollary 2. For I > 1 

(34) gp t = 2gV 2 £ Tl (2*-) + l±t^ (g - 1) |T| 



3=1 



2^q 



where Xj are eigenvalues of the averaging operator T. 

In particular since we know the highest eigenvalue Amax = Q + 1, we immediately 
get (for q > 2) from pijl the asymptotic law for numbers of geodesic paths: gp/ ~ 
as / — > 00 if — q — 1 is not an eigenvalue. If — q — 1 belongs to the spectrum of T, the 
asymptotic is gp; ~ (1 + (—1)') q l ■ For comparison the numbers of all paths behave 
like pi ~ (q + l) 1 . 

4.3. Case q — 1. This is a simple example where we can check (130(1 explicitly. The 
graph T is necessarily a polygon with L > 3 angles. Each long geodesic 7 is character- 
ized by the winding number r > 1 and its direction of movement around the polygon: 
either clockwise or counterclockwise; A(7) = L for all 7. The spectrum of T can be 
found easily: Xj = 2 cos ^ and the trace formula (|5Uj) gives 

L 00 00 

(35) Z T (t) = Y^e 2tcos2 -? = LI (2t) + 2L^I rL {2t) = L ^ I rL {2t) 

j=l r=l r=— 00 

This identity can checked by Poisson resummation. In the limit L — ► 00, t = L 2 r 
(keeping r fixed) we recover a special case of the modular transformation for Jacobi 
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theta function: 

e 4t 



(36) V e"*^ = -t= Y 

j=—oo r=— oo 

4.4. Remark. The original trace formula for regular graph in j2] (rewritten in our 
notations and for special case of trivial character of the fundamental group of T) states 
that for any function g : Z — ► C such that #(n) = —g(n) for all n G Z and 

oo 

(37) ]>>H|<f /2 <oc 

n=l 

the following holds: 

(38) Eto) = l r l-2^/ ^)^^ + E A (t)^ N/ ^(ItI) 

where g(z) = Xl^L-oo g{ n ) z ~ n and numbers ^ are defined by eigenvalues Xj by equation 
A = v ^(z + ^- 1 j. 

Formula ()30|) follows from (|3*Hj) if we choose g(n) = I n (2y/qt) with corresponding 
g(z) = e f ^ z+z \ Integrals representing the contribution of contractible paths convert 
into one another with the change of variables s = Jq (z + z~ Y ). 

4.5. Remark. The trace formula (|30jl is actually valid for graphs with multiple edges 
and loops. The condition of absence of multiple edges and loops was chosen to simplify 
the combinatorial constructions. 
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